Hyperfine interaction induced decoherence of electron spins in quantum dots 
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We investigate in detail, using both analytical and numerical tools, the decoherence of electron 
spins in quantum dots (QDs) coupled to a bath of nuclear spins in magnetic fields or with various ini- 
tial bath polarizations, focusing on the longitudinal relaxation in low and moderate field/polarization 
regimes. An increase of the initial polarization of nuclear spin bath has the same effect on the de- 
coherence process as an increase of the external magnetic field, namely, the decoherence dynamics 
changes from smooth decay to damped oscillations. This change can be observed experimentally 
for a single QD and for a double-QD setup. Our results indicate that substantial increase of the 
decoherence time requires very large bath polarizations, and the use of other methods (dynami- 
cal decoupling or control of the nuclear spins distribution) may be more practical for suppressing 
decoherence of QD-based qubits. 

PACS numbers: 75.10.Jm, 03.65.Yz, 03.67.-a, 02.60.Cb 



I. INTRODUCTION 

Quantum dots (QDs) are very promising candidates 
for future implementation of quantum computations: an 
electron spin in a QD is a natural two-state quantum 
system, which can be efficiently manipulated by the ex- 
ternal magnetic fields and gate voltagesi*^. Also, the 
QD-based architectures are potentially scalable, and rely 
on well-developed semiconductor technologyi. However, 
due to interaction with the environment, the electron spin 
loses coherence very quickly, on a time scale of order of 
nanoseconds for typical GaAs QDs^. It is vitally impor- 
tant for realization of QD-based quantum computing to 
understand the decoherence dynamics in detail, in or- 
der to find practical ways of decoherence suppression. 
Moreover, decoherence of open systems is of fundamen- 
tal interest for understanding of the quantum phenomena 
taking place in mesoscopic systems^i^i^, and therefore at- 
tracts much attention from scientists working in the areas 
of nanoscience, spintronics, and quantum control. 

Among different sources of decoherence relevant for an 
electron spin in a QD, the decoherence by the bath of nu- 
clear spins (spin bath) is dominant for magnetic fields less 
than a few Tesla, and experimentally relevant tempera- 
tures of tens or hundreds of milliKelvin. Much research 
has been focused on the case of a large external magnetic 
field or large bath polarizations, where the perturbation 
theory allows an extensive analysis^i^i^ii^. But the inter- 
esting and experimentally relevant regime of moderate 
magnetic fields and/or moderate bath polarization has 
received much less attention. 

Below, we study in detail the influence of moderate 
magnetic fields and bath polarization on the longitudinal 
decoherence of a single spin in a quantum dot, and two 
spins located in neighboring QDs, where the perturba- 
tion theory is not applicable. In this regime, as the mag- 
netic field or the bath polarization increase, the dynamics 
of the electron spins undergoes a transition from sim- 



ple overdamped decay to underdamped oscillations, and 
these oscillations can be detected using existing experi- 
mental schemes. In this transition, the increase of bath 
polarization affects the decoherence dynamics in exactly 
the same manner as the increase of the magnetic field. 
Our results also show that suppression of decoherence re- 
quires very large bath polarizations, suggesting that the 
use of such methods as dynamical decoupling^ i i ^ d ^ ^ or 
control of the nuclear spins distribution^^, may be more 
practical. 

Theoretical description of the spin-bath decoherence is 
a very complex problemi^, where strong correlation and 
essentially non-Markovian bath dynamics play an impor- 
tant role^ d^d^ . In contrast with the well studied boson- 
bath decoherence^!^, decoherence by a spin bath has not 
yet been understood in detail, especially in moderate ex- 
ternal magnetic fields. 



II. SINGLE ELECTRON SPIN IN A QD 

An electron spin in a QD interacting with a bath of 
nuclear spins is described by the Hamiltonian which in- 
cludes the Zeeman energy of the electron spin in the 
external magnetic field Bq and the contact hyperfine 
couplingiLi^: H ^ HoS, + Ef=i^fcS • Ifc, where S 
is the operator of the electron spin, 1^ is the opera- 
tor of the fc-th bath spin (fc = 1,2,..., N), and Ak = 
(87r/3)5|/iBg„/i„u(xfc) is the contact hyperfine coupling 
which is determined by the electron density u(xfc) at the 
site Xfc of the fc-th nuclear spin and by the Lande factors 
of the electron g* and of the nuclei g„. The terms omit- 
ted in the above Hamiltonian, such as the Zeeman energy 
of the nuclear spins, the anisotropic part of the hyperfine 
coupling, etc., are small, and can be safely neglected at 
the nanosecond timescale, which is considered here. 

In spite of apparent simplicity of the Hamiltonian, 
it is very difficult to determine the dynamics of the 
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central spin S(t). Previously, the quasi-static approxi- 
mation (QSA) for the nuclear spin bath, which treats 
J2k ^felfc a.s a constant, has often been invoke d 
Analytical calculations beyond the QSA have been car- 
ried out only for large magnetic field and/or large bath 
polarizations^iSiiSiiiS, where the perturbative approach is 
valid. Although many qualitative arguments support 
QSA, its validity has not yet been checked in detail. Be- 
low, along with explicit analytical solutions, we provide 
direct verification of our analytics by employing the exact 
numerical simulations for medium-size {N = 20) baths, 
and approximate simulations for large {N — 2000) baths. 
The agreement between all three methods ensures us that 
our findings do not depend on the approximations in- 
volved. 

The analytical calculations within the QSA can be 
performed in different ways which all give the same re- 
sults. The conceptually simplest way is to assume a uni- 
form electron density in the QD, so that all A^. are the 
same, Ak = A, and the Hamiltonian can be written as 
H = HqS^ + AS ■ I, where I = Ifc is the total nuclear 
spin of the bath. This Hamiltonian can be analyzed ex- 
actly, since and Iz + Sz are the integrals of motion. 

First, we consider an unpolarized bath, with the den- 
sity matrix pb{0) = (l/2^)li (g) • • • (g) Iat, where 1^ is 
the unity matrix for k-th bath spin. This density matrix 
can be re- written in the eigenbasis of I^, Iz as Pb{0) = 

Efio Em=-i Pil, M)\I, M) {I, M\ where M = As- 
suming that the initial state of the electron is | t)) the 
quantum-mechanical average of the z-component of the 
electron spin at time t is 



az{t) ^ 2{Sz{t)) =Y.ni,M)?-±^^ (1) 

I,M 



where C = A^ {I - M){I + M + I), B ^ Hq + A{M + 
1/2), and = -I- C^. Due to the symmetry of the 
problem and the initial condition, axit) = cry{t) = 0, 
so we omit these components. The distribution function 
P(/, M) can be calculated^^ and for large and / it 
is approximated by a Gaussian distribution P{I, M) « 
{I/D^/2t:D )e--f'/2^ where D = Replacing the 

summation by integration in Eq. ([1]), we find 



az{t) = l-2W{\D-t), 
and the function VF(A, D] t) has the form 



(2) 



W = l^-^-^^ ^ COS At 



2 2A3 



H — - ^ — + 2$(^=) 



(3) 



where $(2;) is the error function. In this equation, we 
took A = 1, and introduced the notation A = Hq/A; 
this corresponds to normalized energy and time scales, 
so that the time t is measured in the units of 1/A. The 
dynamics of <Tz{t) for several values of A are shown in 
Figs. [Hand [2 
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FIG. 1: (Color online) Electron spin decoherence in various 
magnetic fields [(a), (b), (c)] and with polarized initial nuclear 
spin baths [(d), (e), (f)] for A'^ = 20. The red dashed curves 
denote the exact numerical simulations results for A*' = 20, 
and the blue solid curves correspond to the analytical results. 
The numerical results agree well with the analytical predic- 
tions, and the underdamped oscillations appear once A or k 
is larger than y/N. 



Before discussing the above results, we consider next a 
polarized bath, assuming that its initial state is described 
by a density matrix pb{0) = (l/Z/s) exp{—(3M), where (3 
is the inverse spin temperature, and — [2 cosh(/3/2)]^ 
is statistical sum, and the initial bath polarization is p = 
(M) I {N/2) — — tanh(/3/2). For large spin temperatures, 
the polarization is small, p ~ — /3/2, and we approximate 
P{I,M) « (//i?V2^)e-^'5'/2g-/V2Dg-/3Af^ go ^^^^ 

az{t) = l-2W{K,D;t), (4) 

where k = Df], and the function W{k, D; t) is defined by 
Eq. ([3]). The dynamics of (Jz{t) for several values of /3 is 
shown in Fig. [1] 

The functional form of Eqs. ^ and ^ is identical, up 
to replacing n by A, so the small nonzero bath polariza- 
tion p affects the central spin in exactly the same way as 
the external field of the magnitude Hq — —pAN/2, equal 
to the average Overhauser field exerted on the central 
spin by the nuclear bath. Indeed, the noticeable average 
magnetization (M) = pN/2 of the polarized bath leads to 
a noticeable average Overhauser field, but the variation 
of the magnetization ((AAf)^) = (1 — p'^)N/A changes 
very little at small p, and, correspondingly, the spread in 
the Overhauser fields is almost unchanged. 

From Eqs. ^ and ([3]), we see that (Tz{t) always de- 
cays with characteristic time T2 = \/8/ [NA^). E.g., 
for Ho — (or p = 0), we reproduce a well-known re- 
sult az{t) = (1/3) + (2/3)(l - Dt^)exp{-De/2), i.e. 
az{t) first falls to the value « 0.036, then increases and 
saturates at 1/3. However, for Hq ^ {or p ^ 0), 
due to the oscillatory terms in Eq. ([3]), az{t) can ex- 
hibit oscillations, provided that A (or k) is compara- 
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FIG. 2: (Color online) The same as Fig. [T]but for iV = 2000 
nuclear bath spins in zero magnetic field (a) and in A = 2\/7V 
(b) with initially unpolarized bath. 

ble or larger than \/lJ ^ \/N . This transition, from 
smooth decay at small field/polarization, to the oscilla- 
tions at larger field/polarization, is similar to the well- 
known transition in the dynamics of a damped oscillator: 
the evolution of the central spin is overdamped (or un- 
derdamped) depending on whether the decay time T| is 
larger (or smaller) than the "bare" oscillation frequency 
determined by A or k. Note that the crossover value for 
magnetic field (polarization) is of order of \/N , beyond 
the range of applicability of the perturbation theoryiii^. 

For a typical GaAs QD with the electron delocalized 
over N ^ l{f nuclear spins, A ^ Aq/N 10 ^ /leV 
(where Aq ~ 0.1 meV is the hyperfine coupling for an 
electron localized on a single nucleusSi) . The correspond- 
ing decoherence time Tj* ~ 10 ns (taking into account the 
/ = 3/2 spins of ^^Ga, ^^Ga and ^^As), as confirmed by 
recent experimentsiii^. To observe oscillations for a sin- 
gle electron spin, a very modest external field of order of 
3 mT, or polarization of order of 0.5% is needed. 

It is noteworthy that the decoherence time remains 
practically constant in the course of transition from 
smooth to oscillatory decay. The decoherence time is de- 
termined by the spread in the Overhauser fields, which is 
proportional to the variation of the bath magnetization 
((AM)2) = (1 -p2)7v/4, and it is little affected by smaU 
bath polarizations p^^. For a moderate bath polariza- 
tion, the equivalence between the external field and the 
nonzero bath polarization disappears, and the analytical 
expression for az(t) becomes rather complex, involving 
hypergeometric functions 2^"i(ij b, c, z) of complex argu- 
ment, so we do not present it here. In this regime, the 
decoherence time does not increase much, as shown in 
Fig. WLe) and (f). 

In order to substantially increase the decoherence 
time, an extremely large bath polarization is required, 
as suggested by the results from the perturbation 
methodii^iSiiiS. Such a large bath polarization is beyond 
the scope of our paper. Currently, strongly polarized 
baths are difficult to achieve experimentally, and such 
methods as narrowing the nuclear spin distribution^^ , dy- 
namical decouplings^ and spin echo techniques^'^^, may 
be more practical for suppression of decoherence. 

The quasistatic bath approximation is far from reality. 
E.g., for Ho = and p = 0, the QSA predicts saturation 
of <Jz(t) at 1/3 for t — > oo, which is just an artifact of 
the approximation: the detailed analysis shows that cr^ 



slowly decays (as l/\nt) to zero^ i^^i^" . However, we ex- 
pect QSA to be valid at times of order of j when the 
bath's internal dynamics is not yet important. 

For verification, we perform exact numerical simula- 
tion with medium-size baths of = 20 spins. A real 
QD is approximated by taking nuclear spins located at 
the sites of the 4x5 piece of a square lattice, with the 
lattice constant a = 1. Assuming a parabolic confining 
potential, the electron density is approximated as 2-D 
Gaussian with the widths of Wx — Wy = 1.5, and with 
the center shifted by dx — 0.1 and dy = 0.29 along the 
X- and y-axis respectively. For comparison with the an- 
alytical results above, the effective coupling constant is 

1/2 

defined as A= (^^, ^/N^ ■ The decoherence dynam- 
ics is simulated by directly solving the time dependent 
Schrodinger equation for the wave function of the full 
many-spin system (central spin plus the bath), using the 
Chebyshev polynomial expansion of the evolution opera- 
tor as described in Ref. 

Fig. [TJa), (b), and (c) show analytical [obtained from 
Eq. (12)] and numerical results for (Tz{t), for different mag- 
netic fields. The panels (d), (e), and (f) present az{t) for 
polarized baths with different initial polarizations. For 
the analytical curve in panel (d) the small-polarization 
formula Eq. ^ has been used, while the analytical curves 
in panels (e) and (f) have been calculated from the for- 
mula for moderate initial polarization. The agreement of 
the analytics and the exact numerics is good. The dif- 
ference is caused only by the modest number A'' = 20 of 
the bath spins used for numerical simulations, but exact 
simulations even for A^ = 50 are beyond the capabilities 
of modern computers (since the computation time and 
memory grow exponentially with A^) . 

To study large baths with A^ = 2000, we use the co- 
herent state P-representation described in Ref. [l3. Al- 
though approximate, this numerical approach demon- 
strates excellent accuracy. Figure [5] presents the results 
for az (t) obtained analytically, from Eq. ^ , and numer- 
ically, from P-representation simulations with N ~ 2000 
spins. For A^ = 2000, the agreement between the ana- 
lytics and the numerics is very good. Overall, the data 
presented in Figs. [T] and [2] justify the use of QSA, so that 
we expect our predictions to be valid for real QDs with 
10^ nuclear spins. 

III. TWO ELECTRON SPINS IN DOUBLE QD 

Measurement of the Rabi oscillations for a single elec- 
tron spin in a QD has not yet been achicvcd^^. Recent 
experiments use two electron spins in two neighboring 
QDsiii^. The spins are prepared in the singlet state, 
and the measured quantity is the probability Psit) to 
stay in the singlet state after time t. The oscillations 
described above can be also detected in this double-dot 
setup. When the coupling between the two electron spins 
is negligible, the Hamiltonian of the double-QD system is 
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FIG. 3: (Color online) Time evolution of the probability of the 
singlet state of the two electrons in a double QD by applying 
(a) parallel magnetic fields Ai = A2 = A and (b) antiparallel 
magnetic fields Ai = -A2 = 10^ with Ni= N2= N = 10^. 
Underdamped Rabi oscillations appear once the difference in 
the magnetic fields applied onto each QD is large enough. 



where indices 1 and 2 denote the quantities describing 
left and right QD, respectively (e.g., Hqi and iJo2 are 
the magnetic fields acting on the left and right spins, 
respectively). Note that the evolution of QD 1 and 2 is 
not independent because of the initially entangled singlet 
state, although the Hamiltonian is separable. Using the 
quasistatic approximation, we get 

Ps = lu{XuDi)U{\2,D2) + ^W{Xi,Di)W{\2,D2) 
+ F{X,,D,)F{X2,D2) + G{\i,D,)G{X2,D2) (5) 

where the function W{X,D) is defined by Eq. ([3]), 
andF(A,7:») = 1/2 (^1 + [cos At - (Dt/X) sin Xt]e-^*^ 
G{X,D) ^ 1 - F(A,L>) - W{X,D), U{X,D) = 
il/X'^)[DXtcosXt+ {X^ - D)siiiXt]e-'^''/'^. FigurelSta) 
illustrates dynamics of Ps{t) for unpolarized baths, in the 
case of uniform magnetic field Hf) = Hqi = Hq2- Ps{t) 
decays in the beginning, and saturates at non-zero value 
at long times. The saturation value is 1/3 for zero field, 
and increases with the magnetic field, reaching 1/2 for 
strong fields-'^^. 

However, if the difference between Hqi and Hq2 is com- 
parable or larger than Ay/Ni_2 (which corresponds to few 
milliTcsla for realistic GaAs QD), the probability Ps{t) 
exhibit oscillations [sec Fig.[3^b)], analogous to the oscil- 
lations of (Jz{t) in the single-QD case above. In experi- 
ments, a non-uniform magnetic field can be created e.g. 
by micromagneta^. Another opportunity is using non- 
uniformly polarized nuclear spin baths in the left and 



right QDs. In analogy to the single-QD calculations, it 
can be shown that Psit) in the case of nonzero initial 
polarization is still given by Eqs. ([5]), with replacement 
of Ai^2 by ^1^2- The difference in polarization should be 
about 0.5% for the oscillations to appear. 



IV. CONCLUSION 

In summary, we study in detail the influence of mag- 
netic fields and bath polarization on the decoherence of 
a single spin in a quantum dot, and two spins located 
in neighboring QDs. We focus on the regime of mod- 
erate fields and polarizations, where the perturbation 
theory is not yet applicable, using both analytical tools 
(the quasi-static bath approximation) and the numerical 
simulations (exact, for medium-size baths with 20 spins, 
and approximate, for large baths with 2000 spins). The 
agreement between all three approaches is good, so we 
believe our results are applicable to real QDs with mil- 
lions of nuclear spins. The nonzero bath polarization and 
the external magnetic field influence the decoherence dy- 
namics in exactly the same way, and lead to a transition 
from smooth decay to oscillations once the field (polariza- 
tion) exceeds a certain crossover value. This transition 
can be observed in experiments with a single QD, and 
with two quantum dots. Our results show that substan- 
tial increase of the decoherence time requires extremely 
large bath polarizations, so that such methods as dy- 
namical decouplingiiiiii^ or control of the nuclear spins 
distribution^^ may be more practical for controlling de- 
coherence of QD-based qubits. 
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